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ABSTRACT 

In recent years, a large increase in the number of computer 
controlled systems and processes has dictated the use of devices with 
inherently long time delays. The result has been unstable control 
systems. This thesis develops two methods of feedforward compensation 
which cope with this problem. The Fade' method of polynomial 
approximation is used to formulate a number of transfer functions which 
can be used in the two feedforward compensation schemes. To choose 
the best transfer function for a particular problem, a minimum integrsl 
square error analysis technique is developed for use in the frequency 
domain. Computer simulations are used to verify the effectiveness of 
the two compensation schemes. 
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I. INTRODUCTION 



In any control system, with lumped Linear elements, the output is 
theoretically characterized by a response which starts immediately after 
the application of a transient input. It is well known that this is not 
really true in many cases, since a number of devices may have a 
transportation lag or other delay associated with them. If the delay is 
small, it can usually be ignored. The system is then described by a set 
of linear differential equations and the classical methods of analysis 
can be used to study the behavior of the system. In recent years, tlie 
large increase in the number of com,puter controlled systems and 
processes has dictated the use of devices with inherently long time 
delays. Examples of such devices would be a digital computer controlled 
system, where the computation time could be several seconds, or an 
automated superfractionator, which could have a time delay as long as 
several hours. A typical system, v/ith unity feedback, is shown in Fig. 1. 

When the effects of time delay cannot be ignored, the primary 
problem is one of system instability. This is easily seen by looking at 
the Bode diagram shown in Fig, 2, This diagram represents the open 
loop transfer function for the system in Fig, 1 with 
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PLANT WITH TIME DELAY 
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FIGURE 2. 

Bode Diagram for Equation 1. 
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The time delay, T, was set equal to 1 second. The upper phase angle 
curve represents the system without the time delay. It can be seen 
that the effect of the delay is to increase the phase lag as the frequency 
increases. The amount of the increase is equal to the product of the 
frequency and the time delay. In this particular example, a very stable 
phase margin of 58.5° has been reduced to 24°. If the time delay is 
increased to 2 seconds, the system would be unstable. It can be easily 
concluded that even a modest delay can present formidable stability 
problems. 

This paper presents a method of compensating for time delay in 
closed loop systems, such as the one shown in Fig. 1. The technique 
involves feeding a conditioned signal forward, around the element 
associated with the time delay. The feedforward path may require either 
active or passive elements to condition the signal. The circuit require- 
ments will be developed. This method is an offspring of one referred to 
1,2 

by P. S. Buckley as the Smith Linear Predictor. However, Buckley's 
use of the predictor was not to stabilize a closed loop system. He was 
concerned with the effects of an external disturbance (noise) being 
introduced into the loop between the controller and the plant. He 
assumed the system to be stable before the disturbance v^as applied. The 
predictor was only a method for reducing the effect of the disturbance. 

Thfe Smith Linear Predictor, as presented in Buckley's paper, is shown 
in Fig. 3. 

It can be observed that the circuit requirements for the signal 
conditioner in the feedforward path cannot be synthesized with lumped. 
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linear, passive network elements. A suboptimal network will therefore 
have to bo used in the feedforward path. The Fade' method of polynomial 
approximation will be derived and used to find an acceptable feedforward 
network. It will be found that a number of acceptable approximations 
can be formed by the Fade' method to meet the network requirements. 

The problem will be to choose the best approximation and to find the 
conditions for which it is the best. In this thesis, a frequency domain 
technique is devised for choosing the best approximation in a minimum inters 
gral square error sense. Experimental examples are then given to show 
the effectiveness of the approximations and the compensation method. 
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II. FEEDFORWARD COMPENSATION OF 
SYSTEMS WITH TIME DELAY 



Methods^Coau peji sati on . Given a system such as the one shown in 
Fig. 1, the solution to the stability problem involves the removal of 
the e sT denominator of the closed loop transfer function. 

Two methods of doing this immediately come to mind and ore shown in 
Fig. 4. The transfer function for system 1 (Fig. 4a) can be expressed 



CTs) 

R(S) 



S((s) ^ ^ 




(3) 



It can be seen that if G^(s) is set equal to 02(1 •-e”^'^) the transfer 
function reduces to 



C(S) 



(^) (s) 



ST 



(4) 



I f 6f (s) Gz (s) 

It IS recognized that this equation represents nothing more than the 



transfer function of a closed loop system (without a delay element in 
the loop) cascaded with the time delay. 

In system 2 (Fig. 4b), if G 2 (s) is again set equal to G 2 (l-e“^'^), 
the transfer function can be immediately written as • 

C(s) _ 6,(5} G pJs) 

f?(st ~ I -I- &,(s)6^(s) 

It can be seen that both techniques meet the objective of removing the 
“ST 

e term from the denominator of the closed loop transfer function. 
Both systems use a standard network, v^hich has a transfer function 
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l-e , cascaded with a 02(2) network. This means that the time 
delay network needs to be synthesized only once, normalized v/ith 
respect to the time delay, for all problems. The only requirement is an 
impedance match between the time delay network, 02(2) network, and 
any other sections in the feedfonvard path. There is a subtle difference 
in the two methods of compensation however. This difference involves 
the nature of the output. The two systems can be reduced to the 
equivalent block diagram form shown in Fig. 5 . The output of the first 
system is still delayed by a factor T while the output of the second 
system starts to respond immediately after the application of the input 
signal. First impressions would indicate that the second system is 
better than the first since the time delay has been completely eliminated. 
However, it remains to be seen from experimental data that this is the 
case. It is possible that noise attenuation characteristics, speed of 
response, or other desirable parameter might require' 
use of the first system, A practical consideration might also dictate 
the use of the first system. If, for example, the time delayed plant is 
a chemical process, it would be impossible to "mix" the chemical 
output with an, electrical feedforward signal. On the other hand, the 
electrical feedforv/ard signal could be added to the electrical feedback 
signal coming from some output sensor, as required in the first system. 
The Approximation Requirement. The first step in implementing either 
of these tv/o systems is to synthesize a network which generates the 
transfer function l-e , which is identical to the transfer function 
required for the Smith Linear Predictor, Since it is a transcendental 
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function, some form of approximation must be used if the transfer 
function is to be realized by lumped passive components. Obviously 
distributed parameter delay lines can be used but if T is large, as can 
be expected in a process control, these devices are not practical. A 
set of linear filter transfer functions can be derived using any of the 
standard approximation methods. The background, procedure, and 
accuracy behind the Fade' method of polynomial approximation is 
contained in the next chapter. Approximations based on this method are 
used in this thesis for the feedfonvard network required in the tv/o 
compensation schemes of Fig. 4. 

Summary . The results of this chapter can be summarized as follows: 

1. A closed loop control system with time delay can be stabilized 

if the proper signal can be fed forward around the time delay element. 

2. The filter required to generate the feedforward signal has the 
transfer function 

G(s) = 1-e"^'^ 

3. Approximation techniques must be used to synthesize the 
required transfer function with lumped, linear components. 
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Ill, DESIGN OF OPTIMAL FEEDFORWARD FILTERS 



The Fade' Approximation . A required step in the synthesis of any 

function, in lumped, linear, passive, time invariant form, is to insure 

that the function is rational. The Fade' approximation is a method of 

representing a polynomial or series of polynomials as a rational 
3 4 

function. ' Assume that some function, G(s) , is to be approximated 
by a rational fraction of an n-th order polynomial in s, divided by an 
m“th order polynomial in s. G(s) could, for example, be the Taylor 
series expansion of the desired transfer function Thus G(s) 

can be represented, in general rational form, as 

Q '\r ^ S ^ S 

rrnii ■ «i ■■ 1 1 iiian manirni w ■■■! ( ^ ] 

bo-i- b,S-J- ‘ ‘ ' -i-bm S’’'' 

This expression may bo expanded by cross multiplication to give 




(bot'b, 5 '/" ' 'j( O q -t BiB t " ' ) 

' ( 7 ) 

where G^,G^, etc. are the coefficients of the Taylor series expansion 
of G(s) around the origin in the s -plane. By equating equal powers of 
s, a set of simultaneous algebraic equations are obtained. For a given 
number of terms in the Taylor series and for a particular number of terms 
in the rational fraction, there is a finite number of equations v/hich are 
solved simultaneously to obtain the a's and b's. 
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From the s® terms Aq = 



s 



1 




s 



2 



^3 “ ^3*^0 ^2*^1 ■^' ^1*^2 ^0*^3 



s 



3 



s 



n 




( 8 ) 



The rule for finding the s^ equation should be obvious. It is noted that 

it requires m+n+1 terms of the Taylor series to find n+1 (or n) coefficients 

for the numerator and m (or rn+1) coefficients for the denominator of the 

rational fraction, and that rn+n simultaneous equations are used from (8). 

This still allows one of the coefficients, say bg, to be picked 

arbitrarily and the only difference in the approximations will be a 

constant multiplier for different values of b^. It is realized that when 

the rational fraction is divided out, in ascending powers of s, the first 

m+n+1 terms of the Taylor series are reproduced exactly. Therefore, 

the "maximum error" incurred in using the Fade' approximation is no 

more than the remainder of the Taylor series, taken after m+n+1 terms. 

Nomenclature. The follov/ing nomenclature will be used in the remainder 

of this thesis to refer to a particular approximation or transfer function. 

-sT 

Ej^j^ refers to the Fade' approximation, of the function e , characterized 

by an n-th order polynomial divided by an rn-th order polynomial in s. 

_sT 

But e is not the transfer function required in the feedforv/ard path of 

the systems in Fig. 4. Therefore, will refer to the approximcition 

~sT 

of the desired transfer function, 1-e ", derived from the Fade' 
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approximation of e 



of the same order. For example, if 
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Derivation of the Desired Filter . The Taylor series expansion of e~^, 
about the value x = 0, is the well known power series 



e 



“X 



I-X+ 



Zi 
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( 11 ) 



This is the series that is used to derive rational approximations for the 

-sT -X 

function 1-e , with sT = x. Fade' approximations, for the function e 

are tabulated in many sources.^ Apj^endix 1(a) gives such a list of 

approximations. Appendix 1(b) gives the derived filter transfer 

functions derived from the appropriate approximation, while Appendix 1(c) 

gives the derived filter transfer functions in factored form, Vv^hich is not 

available in the general literature. The polynomials were factored using 

the POLYRT subroutine on the CDC 1604 digital computer and rounded 

off to three decimal places. While the factored form is never used 

directly in any calculations in this paper, it does give a little insight 

as to the location of the poles and zeros of each of the filter transfer 

functions. The poles and zeros are plotted in Fig. 6. There appears to 

be a pattern to the poles and zeros as the order of s in either the 

numerator or denominator is increased. This pattern is not known but 

if it were, it could be used to formulate a general nile for finding the 
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filter transfer function of any desired order. However, the solution 
of (8) is quite straight forward in solving for any particular so that 
no effort was expended to determine an expression for this pole zero 
pattern . 

There is no guarantee that all of the transfer functions found in 
Appendix 1(b) are realizable with passive elements. The properties of 
a transfer function, made up of lumped, linear, passive elements are 
v/ell known but are reiterated here.^ 

1. Poles are in the left half of the s plane or, for the lossless 
case, on the imaginary axis and simple. Complex poles must be 
conjugate pairs. 

2. Zeros may have any s plane location but complex zeros must 
occur in conjugate pairs. A minimum phase transfer function has 
all of its zeros in the left half of the s plane. 

3. Poles may not be located at the origin or at infinity. 

4. The degree of the numerator may equal but not exceed the 
degree of the denominator. 

Item 4 above rules out all transfer functions where m is less than n. 
These transfer functions will no longer be considered in this thesis. 
Referring back to Fig. 6, it is seen that not every Pade' derived transfer 
function can be synthesized with lumped, linear, passive elements, 
even though m is greater than n. As the order of s in the denominator 
is increased, holding the order of s in the numerator constant, the 
motion of the poles is such that poles begin to appear in the right half 
of the s plane. This violates rule 1 above and renders these transfer 
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functions unacceptable as an approximation to the desired filter 
transfer function, 1-e " , It is also noted that for every transfer 
function whei'e m equals n, an amplifier is required, for D.C. gain, in 
addition to the passive elements for synthesis. 

There is no indication as to which of the filter transfer functions 
is "the best." Intuitively, G 77 should be the best approximation to the 
ideal transfer function since its derivation required more terms of the 
povv^er series than any of the others in the list. On the other hand, 
requires more effort to synthesize since it is more complex than the 
others. It may be that a simple filter will work just as well as a 
complex one in the feedforward path of the systems, in Fig. 4. A 
measurement of which is the best approximation to the ideal filter 
transfer function , l-e"^'^, is attempted below. 

The Minimum Integral Square Eri'or Index. The index, chosen to measure 
which of the filter transfer functions is the best approximation, is 

the Minimum Integral Square Value of the difference between the 

"sT 

responses of the ideal transfer function, 1-e , and the Fade' derived 

filter approximation The integral square error can be calculated in 

either the time or frequency domains. Although the time domain has 
received predominant attention in the literature, several advantages 
uncovered during this investigation urge the use of the frequency domain 
technique. The first, and most important of the advantages, is the 
ability to solve for a frequency spectrum which can be later used to find the 
integral square error for any input. This spectrum, called the weighting 
function, gives the user a certain feeling for the "quality" of the 
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network under consideration. When time domain techniques are used, 



an input signal must immediately be chosen and the result applies for 
that input only. The result that is obtained is only a number and is 
not a vei*y satisfying "answer. " Without other networks (numbers) for 
comparison, the results are next to useless. Even with several networks, 
the user only knows which one has the smallest integral square error, not 
what causes that error or what effect another input might have. This 
advantage will be more apparent by comparison of the results of the 
two techniques. Another advantage of the frequency domain approach 
(or disadvantage of the time domain method) is not having to take the 
inverse LaPlace transformation of the network transfer function. Some 
of the larger Fade' approximations are very difficult and time consuming 
when it comes to obtaining time domain solutions. 

In view of the above arguments, this thesis will primarily use the 
frequency domain approach to selecting the best Fade' approximation 
for use as a filter transfer function in the compensation systems of 
Fig. 4. It is known that the mean square value of a signal is equal to 

7 

the value of its autocorrelation function at t = 0. 



where referred to as the Power Density Spectrum of the 

signal. It derives this name from the fact that the integral in (12) equals 



mean sq. value of v(t) = 





( 12 ) 
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the average power dissipated in a 1 ohm resistor by the signal 
voltage v(t). 0^^(jw) can be interpreted as a lineardensity distribution 
of the mean square value of the signal along the imaginary axis of the 
s plane. 0 obtained from the relationship 







and then sotting s = jw. 

Consider the hypothetical experiment shown in Fig, 7. It can be 
7 

shown that the power density spectrum of the error signal Y(s) can be 
expressed as 



0yy{s) = H(S)M(-Sj 

- H(s}6f-S) 

•f- G(S)&(~S)0^^(5) (14) 

where 

H (s) - ( 15 ) 

and is l/PJJtimes the LaPlace transfonn of the autocon'elation 

12 (s) is l/2iitimes the LaPlace transfonn of the cross 
correlation function of Vi(t) and V 2 (t), etc. If vi(t) equals V 2 _(t), as 
shovv'n by the dotted lines in Fig. 7, then all of the auto and cross 
correlation functions in (14) are equal and will be designated by a 
single v) (s). It is known that for any arbitrary function F 

c£s» vv 
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Therefore, from (14) with s = jw 




The quantity in brackets in (17) will be referred to as W(jw) 




:s(l~COS ujt) 5W ujt 



so that 






2(1 



( 19 ) 



From the Fade' approximations, Gj^^(s) is of the form 

/}f^ 3 't S ^ ^ /?/ S Ao 

S c & <p 0 ! ^ 0Q 

where m and n depend on the particular approximation being 

considered. Or in simpler terms 

(sj i- 
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where the M's refer to the sum of the terms with even pov/ers of s and 
the N's refer to the sum of the terms with odd powers of s. For real 
frequencies (when s = jw) 



:r 



( tif) 




( 22 ) 



and 



S(fUJ) 



a 



^ B 



a 



From (18) and (22), it follows that 



(23) 



cos vjfl ’i-Blv4[siN uy'^ 

Substituting (19), (23) and (24) into the bracketed portion of (17) yields 

W(Ji!i)sa(i~cc‘r> ( 25 ) 

“2 (i-ces (vT} ■tB(i’j)SW 

In shortened form, (17) becomes 




W(jw) is given the name Spectral Weighting Function . Referring to 
Fig. 7, the spectral weighting function does just what its name implies. 
It weights the magnitude of the power density spectrum of some input 
signal v(t) to yield the power density spectrum of the error. All that is 
now required, to determine v/hich of the Fade' derived filter transfer 
functions is the best, is the substitution of the various approximations, 
Gmn/ Irito (25) to determine W(jw). Multiplying the spectral weighting 
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function by the power density spectrum of an appropriate input signal 
yields the power density spectrum of the error for that particular 
transfer and input function. Integrating the power density spectrum 

according to (12) yields the integral square value of the error function. 
The Fade' derived filter transfer function that has the smallest area 
under the spectrum, i.e. the minimum integral square error, is the one 
which is the best approximation to the ideal transfer function, l~e . 
The Spectral Weighting Function . In order to generalize the selection 
process for any input signal, the spectral weighting functions, for each 
of the acceptable Fade' derived transfer functions, should be presented 
in graphical form. This permits the user to know in advance which 
frequencies of the input signal the various transfer functions weight 
the most. The transfer function, that best suppresses the frequencies 
in the input signal, is the one that will have the minimum integral square 
error. It is significant to note that the best approximation depends 
upon the spectral distribution of the input signal as well as the spectral 
weighting function. Due to the large number of Fade’ derived filter 
transfer functions and the complexity of the spectral weighting function 
equation, these curves were computed using a digital computer. A 
program was written to compute and plot W(jw) and is presented in 
Appendix 2(b). The results are shown in Appendix 3. The horizontal 
axis for the graphs in Appendix 3 is plotted as the logj^Q of a nonnalized 
frequency time delay product, i.e. x = log^Q wT. In general, the 
spectral weighting functions for the Fade' derived filter transfer 

functions have a larger weighting effect at higher frequencies and 
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appear to act like high pass filters of cutoff frequency, Wq. With 
the exception of the curves corresponding to the case when m equals n, 
all of the spectral weighting functions settle to a constant amplitude of 
one. The spectral weighting function of a transfer function having m 
equal to n has a constant amplitude oscillation about the value two. 

If the frequency components of an input signal were lower than some Wq, 
where the spectral weighing function is zero, then the integral square 
error would be zero. For such an input signal, the Fade' derived filter 
transfer functions would be an excellent approximation to the ideal 
transfer function. 

Comparison of Filter Transfer Functions . The input signal that is quite 
often used to test the performance of a control system or filter is the 
unit step. The idea is that any signal can be approximated by a sum 
of little steps. Any filter or system that would perform well with a 
single step input should work for a sum of steps. A unit step may not, 
however, be a very good choice for an input signal in the selection 
process described above. Note that the minimum integral square error is 
a nonlinear performance index and the principle of superposition cannot 
be used. The Fade' derived filter transfer function that is selected as 
the best, based on a unit step, cannot be assumed to be the best for 
some other input. The minimum integral square error performance index 
does not, however, rule out the possibility that the best transfer 
function for a step input may also be best for other inputs. It all depends 
on the power density spectrum of the input. For the purpose of this paper, 
a unit step will be used as the input signal when selecting a Fade' 
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approximation to use as the feedforward filter transfer function. This 



choice v/as made simply because the unit step is general and easy to 
work with. In an actual problem, where the input signal is known, the 
power density spectrum of the actual input would be used. Also, the 
unit step makes the examination of experimental results a familiar task. 
The power density spectrum of a unit step is given by the equation 



The digital computer is also used to multiply the spectral weighting 
function by the power density spectrum of a unit step and to numerically 
integrate the resulting area under the power density spectrum of the 
error. The program that was written to perform this operation is given in 
appendix 2(c). The integral square error for each of the acceptable Fade' 
derived transfer functions is tabulated in table 1. The table is broken 
into several bands to show how the frequencies of the spectrum cont- 
ribute to the integral square error. Figs. 8 and 9 are graphical 
representations of the data in table 1. The solid lines in Fig. 8 show the 
variation in the integral square error as one of the indices (m or n) is 
varied. Horizontal movement along one of the solid lines from one 
point to another represents a change in the order of s in the numerator, 
while holding the order of s in the denominator constant. It is observed 
that, as long as n is less than m, an increase in the order of s in the 
numerator results in a decrease in the integral square error. The amount 
of the decrease becomes smaller as n approaches m. When n becomes 
equal to m, there is a substantial increase in the integral square error and 
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(27) 
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Table 1 



Band 

Band 

Band 



Band 

G 10 
G 11 
G 20 
G 21 
G 22 
G 30 
G 31 
G 32 
G 33 
G 40 
G 41 
G 42 
G 43 
G 44 
G 51 
G 52 
G 53 
G 54 
G 55 
G 61 
G 62 
G 63 
G 64 
G 65 
G 66 
G 72 
G 73 
G 74 
G 75 
G 76 
G 77 



Integral Square Errors (Frequency Domain) for 
Fade' Derived Filter Transfer Functions 
Step Input Time Delay 1 Second 



less 


than .01 


rps 


Band 3 


.01 to 1 . rps 




Band 4 


1. to 10. rps 




Band 5 


0 


1 


2 


3 


aprx. 0 


.0165 


. 1868 


.0293 


aprx. 0 


.0003 


.2021 


.0717 


aprK, 0 


.0016 


.1142 


.0286 


aprx. 0 


.0000 


.0748 


.0283 


aprx. 0 


aprx. 0 


.0974 


.0506 


aprx. 0 


aprx. 0 


.0974 


.0206 


aprx. 0 


aprx. 0 


.0449 


.0203 


aprx. 0 


aprx. 0 


.0369 


.0290 


aprx. 0 


aprx. 0 


.0401 


.0605 


aprx. 0 


apnc, 0 


. 1474 


.0206 


apnc . 0 


aprx. 0 


.0347 


.0286 


aprx. 0 


apix. 0 


.0197 


.0296 


aprx. 0 


aprx. 0 


.0137 


.0340 


aprx. 0 


aprx. 0 


.0095 


.0695 


aprx. 0 


aprx. 0 


.0366 


.0280 


aprx. 0 


aprx. 0 


.0138 


.0292 


aprx. 0 


aprx. 0 


.0064 


.0306 


aprx. 0 


aprx. 0 


.0030 


.0344 


aprx. 0 


aprx. 0 


.0012 


.0584 


aprx. 0 


aprx. 0 


.1366 


.0286 


aprx. 0 


aprx, 0 


.0130 


.0206 


apr>:. 0 


aprx. 0 


.0040 


.0284 


aprx. 0 


aprx. 0 


.0012 


.0292 


aprx. 0 


aprx. 0 


.0002 


.0300 


aprx. 0 


aprx. 0 


.0000 


.0486 


aprx. 0 


aprx. 0 


. 0234 


.0284 


aprx. 0 


aprx, 0 


.0036 


.0272 


aprx. 0 


aprx. 0 


.0006 


.0252 


aprx. 0 


aprx. 0 


.0000 


.0242 


aprx. 0 


aprx. 0 


aprx, 0 


. 0254 


apn< . 0 


aprx. 0 


aprjc. 0 


.0410 



10. to 100. rps 
100. to 1000. rps 
over 1000. rps 



4 


5 


M.S.K. 


.0029 


apr>c. 0 


.2355 


.005 7 


aprx, 0 


.2699 


.0029 


aprx , 0 


.1471 


.0029 


apnc, 0 


.1060 


.0057 


aprx , 0 


.1537 


.0029 


aprx, 0 


. 1289 


.0029 


aprx. 0 


.0761 


.0029 


aprx, 0 


.0607 


.0057 


aprx. 0 


.1063 


.0029 


aprx, 0 


.1789 


.0029 


aprx. 0 


.0662 


.0029 


apnc, 0 


.0522 


.0029 


aprx. 0 


.0508 


.0057 


apnc, 0 


.0009 


.0029 


aprx. 0 


.0682 


.0029 


apnc. 0 


.0450 


.0029 


aprx, 0 


.0398 


.0029 


aprx. 0 


.0402 


.005 7 


apix. 0 


.0652 


.0029 


aprx, 0 


. 1681 


.0029 


aprjc . 0 


.0446 


.0029 


aprx, 0 


.0353 


.0029 


apnc , 0 


. 0322 


.0029 


aprx, 0 


.0332 


.0057 


aprx. 0 


. 0544 


.0029 


aprx, 0 


.0545 


.0029 


apnc, 0 


.0338 


.0029 


aprx. 0 


.0287 


.0029 


aprx . 0 


.0271 


.0029 


aprx, 0 


.0203 


.0057 


aprx, 0 


.0467 
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wh©n n becomes larger than m the integral square error goes to infinity. 
Moving vertically from one solid line to another represents a change 
in the order of s in the denominator while holding the order of s in the 
numerator constant. For these changes, the integral square error usually 
decreases as the order of s is increased. The exceptions to this 
occur when there is a pole "close" to the imaginary axis (Fig. 6) such 
as G4Q and Gg]^. The dotted lines in Fig. 8 represent the change in the 
integral square error as both indicics are changed, but the difference 
is held constant. It is readily obseiwed that the amount of improvement 
in the integral square error is exponentially decreasing as the order (and 
complexity) of the transfer function is increasing. Figs. 9 (a and b) 
show the changes in the contribution to the integral square error by the 
various frequency bands. It is observed, from Fig. 9(a), that changing 
the order of s in the denominator, while holding the order of s in the 
numerator constant, has little effect on which band contributes the most 
to^the integral square error. On the other hand. Fig. 9(b) shows 
that changing n while holding m constant has a considerable effect on 
the frequency band contributing the most area. As n is increased, the 
higher frequency bands tend to contribute more. This fact could be used 
to advantage in many systems. A large number of closed loop systems 
tend to act as a low pass filter. If one of the Fade' derived filters is 
used in such a closed loop system, the low pass filtering effect of the 
system should suppress the higher frequencies that are weighted more 
by the high pass feedforward filter. This means that the integral square 
error for a unit step input to such a system would be less than that 
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when n becomes larger than m the integral square error goes to infinity. 
Moving vertically from one solid line to another represents a change 
in the order of s in the denominator while liolding the order of s in the 
numerator constant. For these changes , the integral square error usually 
decreases as the order of s is increased. The exceptions to this 
occur Vv'hen there is a pole "close" to the imaginaiy axis (Fig. 6) such 
as G^q and Ggj^. The dotted lines in Fig. 8 represent the change in the 
integral square error as both indicies are changed, but the difference 
is held constant. It is readily observed that the amount of improvement 
in the integral square error is exponentially decreasing as the order (and 
complexity) of the transfer function is increasing. Figs. 9 (a and b) 
show the changes in the contribution to the integral square error by the 
various frequency bands. It is observed, from Fig. 9(a), that changing 
the order of s in the denominator, while holding the order of s in the 
numerator constant, has little effect on Vi^hich band contributes the most 
to the integral • square error. On the other hand. Fig. 9(b) shows 
that changing n while holding m constant has a considerable effect on 
the frequency band contributing the most area. As n is increased, the 
higher frequency bands tend to contribute more. This fact could be used 
to advantage in many systems. A large number of closed loop systems 
tend to act as a low pass filter. If one of the Fade' derived filters is 
used in such a closed loop system, the low pass filtering effect of the 
system should suppress the higher frequencies that are weighted more 
by the high pass feedforward filter. This means that the integral square 
error for a unit step input to such a system would be less than that 
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calculated above. It would therefore be advantageous to use the filter 
transfer function with the highest frequency weighting in systems that 
act as low pass filters. 

The Integral Square Error Index, (Time Domain). In order to check the 
above results, the integral square error is calculated for a fevz cases 
in the time domain. The same hypothetical experiment, shown in Fig. 7, 
is used. The equation for the error" signal is 




(28) 



Taking the inverse LaPlace Transfonn 
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The integral square error is given by the equation 
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It is easily seen that the first integral reduces to T (the area under a 
pulse of magnitude 1 and duration T) , The second integral is equal to 
the area under the response curve of the Fade' derived transfer function, 
GjTin/ froni time 0 to T. The last integral represents the area under the 
square of the filter response curve for all time. The lower limits can 
be changed to 0 since the response does not start prior to the application 
of the input at time 0. The equation for the integral square ciTor, for a 
unit step input, therefore reduces to 



The digital computer program presented in appendix 2(d) was 

written to compute I from (33). Because of the long computation time, 

only a few of the Fade' derived filter transfer functions were tested by 

this method. The results are given in table 2. Since all of the 

approximations contain decaying exponential terms, the upper limit for 

the third integral in (33) was set at 20 times the longest time constant. 

By comparing tables 1 and 2, it is seen that the results of the two 

methods agree to within the third decimal place. Better agreement can 

be obtained if the tolerance on the numerical integration is made smaller. 

Such action would increase the computation time however. 

“sT 

Other Approximations for 1-e Several other sets of approximations 

were investigated but none performed as well as the Fade'. The defining 
equations for the other approximations investigated are 




(33) 




( 34 ) 
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Table 2 



G 10 
G 11 
G 20 
G 21 
G 22 



band 

B 1 
B 2 
B 3 
B 4 
B 5 

C 1 
C 2 
C 3 
C 4 
C 5 



Integral .Square Errors (Time Domain) for Several 



Pads' Transfer Functions - Step Input 



Time Delay - I Second 



first 


second 


third 




integral 


integral 


integral 


I.S.E. 


1, 


.6321 


.5000 


. 2358 


1. 


.8647 


1.0000 


.2707 


1. 


.8006 


.7490 


. 1477 


1. 


.8818 


.8696 


.1059 


1. 


.9225 


.9992 


.1543 




Table 3 







Integral Square Errors for 
Bj to Bg and to 
Step Input - Time Delay 1 Second 

Band designations are the same as for Table 1 



0 


1 


2 


3 


4 


5 


M.S.E. 


apioc. 0 


. 0166 


.1868 


. 0293 


.0029 


aprx. 0 


.2355 


aprx. 0 


.0054 


. 1289 


.0286 


.0029 


aprx. 0 


. 1662 


aprx. 0 


.0025 


. 1012 


.0284 


.0029 


aprx. 0 


.1353 


aprx. 0 


. 0016 


. 0840 


.0286 


.0029 


aprx. 0 


.1168 


aprx. 0 


.0010 


. 0719 


.0286 


. 0029 


aprx . 0 


.1044 


aprx. 0 


.1248 


. 1808 


.0608 


.0057 


aprx. 0 


.4676 


aprx. 0 


.2779 


.2852 


.0512 


.0057 


aprx. 0 


. 6201 


aprx. 0 


. 2928 


.4262 


. 0633 


.0057 


aprx . 0 


. 7884 


aprx. 0 


.2986 


.5621 


.0589 


.0057 


aprx . 0 


.9253 


aprx. 0 


.3014 


.6621 


.0484 


.0057 


aprx. 0 


1.0176 



41 



and 




These were chosen mainly on the basis that they are relatively easy to 
synthesize. For example, (35) is a lattice network yielding the proper 
time delay. The results of the investigation are given in table 3. The 
format is the same as that for table 1 (i.c. the frequencies in each 
band are the same as for the bands in table 1) , It is noted that the integ- 
ral square error for the approximation becomes smaller as n is 
increased. However, for the same degree of complexity, the Fade' 
approximation has a smaller mean square error. The Cj^ approximation 
does not perform v^^ell at all. The smallest integral square error occurs 
when n = 1 and as n is increased, the ititegral square error increases. The 
smallest integral square error for a approximation is larger than the 
largest integral square error for a Fade' approximation. 

Summary . From the above observations, it appears that the Fade' 
approximation derived from the largest number of terms in the power 

series has the smallest integral square error and is therefore the best 

-sT 

approximation to the ideal transfer function, 1-e . This statement 

is, of course, subject to the condition that the denominator is of higher 
order than the numerator and that there are no poles in the right half of 
the s plane. It is speculated however, that the use of this approximation 
may not be the most practical in many cases. The increase in complexity 
(and therefore in difficulty in synthesis) may far outweigh the advantage 
of a small decrease in the integral square error. As stated before, the 
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selection of how complex a Fade' derived filter transfer function to use 



will depend on the system being compensated, the duration of the time 
delay, the nature of the input signal, etc. The results of this chapter 
can be summarised as follows; 

1. The Fade' approximation is a method for representing a 
transcendental function as a rational fraction of polynomials. It 
is found by the simultaneous solution of a series of algebraic 
equations. 

“sT 

2. Some of the Fade' approximations of the desired function 1-e 
can be realized with lumped, linear passive elements. 

3. The best approach in choosing a filter transfer function is' 
through the frequency domain. In this procedure, a quantity called 

the Spectral Weighting Function is generated from the desired transfer 
function and the transfer function of a filter under consideration. 

The spectral weighting function graphically shovjs i^hich frequencies 
in the input signal contribute most to the integral square error of the 
filter. Idien the spectral weighting function is multiplied by the 
power density spectrum of an input signal, the power density 
spectrum of the error is obtained. The filter with the minimum area 
under the power density spectrum of the error is the one with the 
minimum integral square error and is considered "the best." 

4: The Fade' approximations have the smallest integral square error 

of all of the approximations investigated for the step input. 
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IV. EXPERIMENTAL INVESTiaATIONS 



In order to examine the effects of approximation in the feedforward 

compensation schemes outlined earlier, a hybrid computer model was 

tested with several of the Fade' derived filters in the feedfoiward path. 

The specific objectives of the tests were to ascertain the derogatory 

effects of using an approximation, rather than the ideal transfer function 
"“ST 

l-e , and to determine v/hat effects the order of the approximation 
has on the time domain performance of the feedforvvrard filter. 

The Plant . For simplicity, the plant described in the introduction was 
the one used for the hybrid computer model. The plant is shown again 
in Fig. 10(a). Fig. 10(b) shows the computer model. An EAC model 
TR-20 analog computer was used to simulate the linear portions of the 
plant. A CDC model 160 digital computer was used as the plant time 
delay element. In addition, a portion of the analog computer provided 
signal conditioning for proper A/D and D/A conversion. The A/D 
converter requires all signal inputs to be between 0 and ~10 volts. The 
analog computer setup is shown in Fig. 11. The maximum allowable 
signal at any point in the computer model is chosen to be +8 volts. 

This is a convenient choice of limits since the digital computer 
operates on the octal number system and a direct correlation between 
voltage levels in the analog computer and numbers in the digital 
computer memory results. 

The Digital Time Delay . The basic idea behind the digital time delay 
is not complex. The analog computer output is sampled periodically 
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and the number corresponding to the amplitude of that sample is stored 
in the computer memory. At the end of the time delay period T, the 
number is returned to the analog computer as a voltage signal through 
a hold element. If the sampling rate is fast enough, the voltage 
returned to the analog computer appears to be a continuous time varying 
signal. A flow graph of this operation is shown in Fig. 12. The 
computer program v/ritten to perform this operation is given in appendix 2 (f). 
There is a difficulty in the use of the program listed in appendix 2(f) 
however. As the time delay is increased, the sampling rate is decreased. 
This program is only adequate for for delays of several minutes or less. 

For time delays longer than this, the output of the D/A converter is better 
described as piecev/ise constant rather than continuous. 

Simulation of a Typical Fade' Derived Filter . The most difficult step in 
setting up a computer model of the system is the simulation of the 
feedforward filter, Large gains (larger than 10^ in G 77 ) are 

required in all but the simplest of the Fade' derived filters. The steps 
taken to simulate G^g are given here as an example of the procedure. 




Dividing the numerator and denominator of (36) by (sT)'^ and setting 



equal to the ratio of the output to the input. 




( 37 ) 
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By cross multiplication and rearranging terms 
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( 38 ) 



(38) is in the proper form for the simulation on the analog computer. 

The analog computer setup for 0^3 is shown in Fig. 13. The other Fade' 
derived filter approximations are simulated in the same way. 

The Step Response . It was stated in the introduction that the system 
under consideration v/as stable with a time delay of 1 second but 
unstable if the time delay v/as increased to 2 seconds. The Bode 
diagram for the plant was presented in Fig. 2 to support this. To verify 
this, the feedforward filter v/as disconnected in the computer model and 
the model was subjected to a 1 volt step. The results are shown in 
Fig. 14. The plant is definitely unstable for a time delay of 2 seconds. 

To meet the objectives of the tests, three of the Fade' derived 
filter transfer functions were used as the feedforward network in the 
computer model. The three that were used are G 2i> G 43 / 0 ^ 5 . 

These v/ere chosen because they had the minimum integral square error 
for a given value of m. The time delay was set to 2, 3, and 4 seconds. 
At all three of the time delay settings, the plant without compensation 
would be unstable. The results of the tests, for both compensation 
methods, are shov/n in Figs. 15 (a, b and c) . It is observed that all 
three Fade' derived filter transfer functions stabilize the system for all 
three values of time delay tested. For the longer delays, the more 
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Step Response with G21 Compensator 
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FIGURE 15(b) 

Step Response with Compensator 



53 





FIGURE 15(c) 

Step Response with Gj:^^ Compensator 
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complex filters do a better job with respect to the shape of the response 
curve. The response with the Gg^ filter and with a 4 second time delay- 
looks just like the response of the uncompensated system with no time 
delay, with the exception of the time shift. It is interesting to note that 
the two compensation methods (Fig. 4) do not give the same compensated 
system output. The best example of observing this is the response of 
G 21 with a 4 second delay, Fig. 15(a). The response due to compensa- 
tion method 2 oscillates much more than the response due to compensa- 
tion method 1 and the settling time is longer. On the vdiole, it appears 
that compensation method 1 performs better than compensation method 2. 
This is very fortunate since it was noted earlier that practical 
considerations might require the use of this compensation scheme. 
Sensitivity . No extensive investigation v/as made into the area of 
sensitivity. Only one short experiment was performed to get an idea of 
the tolerance limits on the filter time delay setting. The plant was 
simulated on the hybrid computer with a two second time delay. This 
is an unstable condition. The plant was then compensated, using 
compensation method 1 and the Fade' derived filter G 2 J. The filter time 
delay was vari.ed from 1 to 5 seconds, in half second steps. The 1 
second setting was the shortest delay that could be conveniently set 
on the computer and at 5 seconds, the system v/as at the margin of 
stability. The results of some of the other trials are shov/n in Fig. 16. 

It is observed that the system is stable for a large range of settings of 
filter time delay. The quality of the response deteriorates as the 
difference betv/een the plant and filter time delays increases but the 
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V, CONCLUSION 



It was easily shown, in Chapter 2 of this thesis, that the 
feedforward compensation techniques presented in Fig. 4 would remove 
the phase lag due to a time delay in a closed loop control system and 
v/ould therefore improve the stability of such a system. It was then 
implied that the required feedfonvard network could not be synthesized 
with lumped, linear, passive, time invarient circuit elements. The 
problem was therefore reduced to one of finding a rational approximation 
that would perform as near to the unrealizable ideal circuit as possible. 
Since the function to be synthesized is a common one, a number of 
approximations and approximation techniques were available to choose 
from. The performance index, chosen to rate the various approximations, 
was the minimum integral square error between the step response of the 
ideal and approximation filters. The Fade' method of polynomial 
approximation was finally used since it was simple to apply and it gave 
a wide variation in transfer function complexity. This large variety of 
circuits means that a large range of integral square error v^as available to 
test the effects of approximation in the tv/o compensation schemes. It 
was found, through computer simulations, that approximations of the 
ideal transfer function effectively stabilize systems with time delay. 

Even poor approximations, i.e. filters with large integral square errors, 
were able to effect a degree of stability in such systems. In addition, 
it was found that the stabilizing properties of the two compensation 
schemes were not affected by small changes in the time delay setting 
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system remains stable. This feature is useful in situations where the 
plant time delay is either not known exactly or it changes. 

Summary, The observations made in this chapter can be summarized 
as follows: 

1. The Fade' derived filter transfer functions make effective 
feedforward filters in the compensation schemes outlined in Fig. 4. 

As anticipated, the more complex transfer functions, which are 
better approximations to the ideal transfer function, give a better 
system response. 

2. Compensation method 1 (Fig. 4(a)) gives a better system response 
than compensation method 2. 

3. The ability of the Fade’ derived filter transfer functions to 
stabilize an unstable plant is not vei~y sensitive to changes in the 
filter time delay setting. 
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V. CONCLUSION 



It was easily shown, in Chapter 2 of this thesis, that the 
feedfonvard compensation techniques presented in Fig. 4 would remove 
the phase lag due to a time delay in a closed loop control system and 
would therefore improve the stability of such a system. It was then 
implied that the required feedforv^ard network could not be synthesized 
with lumped, linear, passive, time invarient circuit elements. The 
problem was therefore reduced to one of finding a rational approximation 
that would perfonn as near to the unrealizable ideal circuit as possible. 
Since the function to be synthesized is a common one, a number of 
approximations and approximation techniques v<7ere available to choose 
from. The performance index, chosen to rate the various approximations, 
was the minimum integral square error between the step response of the 
ideal and approximation filters. The Fade' method of polynomial 
approximation v/as finally used since it was simple to apply and it gave 
a wide variation in transfer function complexity. This large variety of 
circuits means that a large range of integral square error was available to 
test the effects of approximation in the two compensation schemes. It 
was found, through computer simulations, that approximations of the 
ideal transfer function effectively stabilize systems with time delay. 

Even poor approximations , i.e. filters with large integral square errors, 
were able to effect a degree of stability in such systems. In addition, 
it was found that the stabilizing properties of the two compensation 
schemes were not affected by small changes in the time delay setting 
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of the approximation feedfonvard filter. In general, it can be concluded 
that both the compensation methods presented in Chapter 2 are excellent 
compensators for time delay and that the required feedforv/ard circuits 
are not very sensitive to approximation effects. 

The concept of a frequency domain analysis of the minimum integral 
square error criterion, as developed in Chapter 3, is as important as 
the idea of feedforward compensation. In most of the literature, the 
time domain analysis is used. Vvdule the time domain technique gives 
satisfactory answers, i.e. gives the integral square error for a problem 
under consideration, the frequency domain analysis gives this 
information and much more. The spectral weighting function, an inter- 
mediate step in the frequency domain analysis, gives the user a 
graphical indication of the frequencies that are contributing to integral 
square error. This means that the user can find a solution for his 
problem based on the frequency spectrum of the input signal and the 
spectral weighting functions of the various solutions under consideration. 
The spectral weighting function can be computed once and used to 
calculate the integral square error for any input. Ihese advantages are 

not available when the time domain analysis is used. 

Topics for Further Study. The ideas used in this thesis can be projected 
into additional areas of study. For example, work can be done in the 
case where the time delay cannot be isolated from the remainder of the 
open loop transfer function of the plant. Also, the idea of a sensitivity 
investigation could receive additional attention. It is interesting to note 
that the ideal feedforward transfer function is very similar to the 
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transfer function of a Zero Order Hold circuit, (l™e~^'^)/s. A number 
of the Fade' derived filters found in appendix 1 could possibly be used 
in conjunction with an integrator to produce a passive approximation of 
a zero order hold circuit. The Fade' polynomials are a promising method 
for finding a rational approximation for any irrational function and work 
can be done in the implementation of such functions. Another area 
where work can be done is to consider phase relationships in the 
frequency domain analysis of the error. Only the amplitude relationship 

(integral square error) v;as considered in this thesis. It would be 
interesting to see what additional information could be obtained from 
such an analysis. 
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6o-36x-t-9X -X 
2 

60h24x+5X 

2 3 4 

360-2toX+72X -12X +X 
2 



E 



360+120X+12X 

6 



30 



31 



2 3 
6+6X+3X +x 



24-6x 



E 



32 



33 



34 



E 



40 



E 



4l 



2 3 

24+i8x+6x +x 
22 

6o-24x+3X 
2 5 

60+36X+9X +x 

2 3 

120-60X+12X -X 
2 3 

120+60X+12X +X 

2 34 

84o-48ox+12X -16x +X 

^ — 3 

84o+36ox+6ox 44x 
24 

2 — 

24+24X-H2X -t-4x +X 

120-24X 

2 3 

120+96X+36X +8X -i-X 



D = 360-120X-1-12X 

k2 2 5 

360+2iiOX-!-72X +12X +X 
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# 



2 5 

E = 840-'^6ox-f6ox -Ux 

45 ‘ 2 5 

840-!-480X il20X -}l6x +X 

2 54 

^ = -'>-^-640X+18QX -20X -!-X 

44 ‘234 

l680•!-8'40X^■l80X +20X -:X 

E = 120 

50 “2545 

120+120X-1-60X -I20X +5X -!-X 
E = 72 O-I 2 OX 

51 2 3 n~T 

720+6oOX+24ox -i 60X +10X -!-X 
2 

E = 252 O- 72 OX+ 6 OX 

52 23 4 ~T 

252 O+I 8 OOX+ 60 OX -! 120 X -!- 15 X +X 

2 3 

E = 672 O- 232 OX+ 36 OX -20X 

53 23 
6720 + 4200 X+I 200 X + 200 X + 20 X -f-X 

2 54 

E = I 512 O- 672 OX- 1 -I 26 OX - 120 X + 5 X 

54 ~ 3“^ 4 5 

13120 + 8 tiOOX+ 2 lOOX + 3 OOX + 25 X +X 

2 345 

E = 30240 -I 512 OX+ 336 OX -420X + 3 OX -X 

55 " ~2 3 

30240 +I5120X+3360X +420X +30X -i-X 

E = 720 

60 2 3 4 T~6 

720+720X+360X +120X -1-3OX +6x +X 

E = 5040-720X 

61 ~2 5 

5040+4320X+I800X -i-400X ^-90X -M2X H-X 

2 

E = 20180-50»K)X+36oX 

62 • ‘ 234 “5 ^ 

20 l 60 il5.120XH-54O0X + 1200 X -}-l 8 ox +i 8 x -t-x 

2 3 

E = 6o4&0-2016oX+2520X -120X 
63 ' 2 3 5 S' 



6o480+40320X+126oOX +24oOX -!-300X +24x +X 



64 



E 

6k 



E 

65 



E 

66 

E 

70 



E 

71 



E 

72 



E 

75 



2 5 

151200 - 6 oi^- 80 X+ 10080 X - 8 !-: 0 X • 1 - 50 X 

^ 

151200-!-90720X-:-25200X -?-llf-200X -<-’l50X +50X hX 

2 5 5 

35264 o-151200X-t-5024ox - 556 OX ■h23.0X - 6 x 

2 5 ' U 5 ^ 

33264o-i-i8i4Uox+4536ox -!-6720X -i-630X +56x +x 

2 3 456 

665280 - 332640 X+ 73600 X - 10080 X -»64px -4^ +x_ 

2 ' 3 456 

66528o:-3526lK)Xi-756oox + 10080 X -i- 8 t! 0 X -iJt 2 X +x 

504o 

2 3 4,567 

50^,0+5o4ox+2520x -t84ox + 210 X 7-I+2X +7X -ix 

4o320-504ox 

- ^ —2 5 4 5 5 Y 

40320-i-35280X+15120X -tl, 200 X H-840X + 126 X +l4x +X 

2 

l8l4i^0-40320X-l•2520X 

- — 2 3 4 5 6 7 

18i440+14i120X+52920X +I 260 OX + 2100 X + 252 X + 2 K +X 

2 3 

6o48oo-i8 i44ox-i-20i6ox -84 ox 

2 3 5 5~7 

6o48oo-f42336ox+i4n2ox +29400X +4200X +420X -i- 28 x +x 



234 

E = 1663200 - 604800 X+ 90720 X - 672 OX + 210 X __ 

74 23 ' T 5 F~ 7 
i 6632004 - 10584 oox+ 317520 X +583oox +7520X + 630 X +35X +x 

2 545 

E = 399168 O-I 66320 OX+ 30240 OX -302!K)X + 168 OX -42X 

75 23 5 S~T 

3991^0-{-23284^ox+635o4ox +10584ox +L 1760 X + 882 X +fv 2 X +x 

2 3456 

E = 8648640 - 399 l 68 ox^ 63 l 6 oox -IOO 80 OX +756ox -336x +7X 

76 234 5 

864864o-f465696ox+ll6424ox +170-iOOX +17640X +II 76 X ■tl^ 9 X +x 

2 3 4 567 

E = 17297280-8648640X+1995840X -277200X -!-25200X -1512X +56x ::X_ 

YY — - ^ 2 5 4 5 

172972804-864864 ox+ 199584 ox + 277200X +25200X +1512X +56X -iX 
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/iPPEKDK 1(B) 



G 

ran 



-sT 



PAGE' /iPPROXBlATION OF 1-e 



1 E 



X = sT 



ran 



G = 0 

GO 

G = X 
01 

2 

G = 2 X-X 

02 2 

2 5 

G = 6 X- 3 X +X 

05 6 

2 5’^ 

G = 24 X- 12 X ^J^X -X 

o 4 24 

G = _L_ 

10 1+X 

G = 2X 

11 

2 

G = 6 x-x 
12 "5+2x 

2 5 

G = 24x-6x +X 
15 24+^"^ 

254 
G = I 2 OX- 56 X + 8 x -X 
l 4 120454 X 

2 

G = 2 X+X 
20 2 

2 + 2 X+X 

2 

G = 6 x-!-X 

21 2 

6 V 4 X+X 
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G 

22 



G 

25 



G 

24 



G 

50 



G 

51 



G 

52 



G 

55 



G 

54 



G 

4o 



G 

4l 



G 

42 



12X 

2 

12+6X+X 

2 5 
6ox-6x +x 

2 

6o+24x^-5X 

2 5 4 

560 X- 60 X + 12 X -X 
2 

560-fl20XHl2X 

2 5 
6 x+ 5 x +x 
■'?— 5 

6+6x+ 5X ^x 

2 5 

24x+6x +x 

2 5 

24-M8x^-6x +X 

2 5 

6ox-f6x -fX 

2 5 

60+56X+9X +x 

5 

120X+2X 

2 5 

120+60X+12X +X 

2 5 ^ 
8tl0X-60X -{-20X -X 
2 5 

84o+56ox+6ox 4-4x 
254 

24X+12X hJiX +X 
2 

24+24X-I-12X 4-4x +X 

254 
120X+56X +8x H-x 
“^^ 5 *+ 

l20 i-96x-f 56X +8x -tx 

2 34 

56 OX+ 6 OX -{12X +X 

2 5 ~^ 

360+2!-!GX4-72X -M2X ^•X 
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G 

45 



2 5 

.0_4DX+6oX -! 20X +)C 

“~2 5"ir 

8!^0-f460X+120X +l6x -i-X 



5 

G = i68ox-!4ox 

44 ~2 5 4 

l63o+84ox-^l8ox -1-2OX +X 



G 

50 



2 545 



I2OX+6OX +20X +^X +X 
' 2 

120-i-120X-i8oX +20X +5X +X 



G 

51 



G 

52 



G 

55 



G 

54 



2545 

720X+240X +6OX -!-10X +X 

Y 5 ITT 

720+600X+24oX + 6OX +10X +X 

2 545 

2520X-I540X -H20X -t-15X +X 

2 's~^r 5 

252O+18OOX460OX +120X +15X -l-X 

2 545 

6720X-i84ox -<-220X +20X 4X 

2 5 4 5 

6720+4-200X4-1200X +200X +20X +X 

2 545 

15120X-t840X -i420X +20X +X 

2 5 

15120+8iK)CK+2lOOX +5OOX +25X +X 



5 5 

G = 5024ox-f84ox -h2X 

55 2 3 4~5 

502^^0+15120X+5560X +420X + 50 X +x 

2 5456 

G = 72 OX+ 360 X + 120 X + 50 X +6x +X 

60 ~ 5 4 5 S’ 

72 O+ 72 OX+ 56 OX + 120 X + 50 X 46 X +X 

2 3456 

G = 5o4ox+i8oox -!-48ox + 90 X +12X +x 

61 25 4 . . 5 ~T 

5040-f4520X+l800X -fj| 80 x + 90 X + 12 X 4 -X 

2 5456 

G = 20160X-i- 5040X 4-1200X 4-l80X 4-l8X 4X 

62 2 ^ 4 5 " () 

20160+15120X4-5Ji 00X +1200X -1-I8OX +l8x +x 
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2 3 4^6 

G = 60480X-il0080X +2520X +300X + 2 Ux +X 

^ 

60480^-'^0520X-i•12600X •;-2400X +500X + 2 Ux -iX 

2 5 U 5 6 

G = I312OOX+3512OX +3040X 4t20X H-30X +X 

64 2 3 4 ” 5 Z ~ 

I5I2OO+9072OX+2520OX +4-200X ^J350X +30X +X 

2 5456 

G = 332640X-H5120X H-IOOSOX -f420X H-42X +X 

65 2 ^ 3 4 5 6" 

332640+1814110X+45360X +6720X +63OX +36X +X 



G = 

66 



G 

70 



G 

71 



G 

72 



G = 

73 



G 

7^- 



G 

75 



G = 

76 



66328OX+20160X -f84x 

_ _ ^ ^ g. 

665280+332640X+75600X +IO08OX h84ox +42X +X 
2 3 4567 

304ox+2320X 4840X +210X +42X h-7X -<X 
2 ~5 4 5 6~T 

504o+5o4ox+2520x +84ox +210X H-42X +7X +X 

2 54567 

4032OX+I512OX •f4200X 484 ox -^-126x -H4x +X 

2 3 4 5 r~J 

40320+35280X+I5120X +4200X +840X +126X +l4x +x 

2 3 4567 

18i 440X+504Q0X -Q26OOX +2100X -i252X +21X +X 

2 5 4 5 6 7 

18i44o+14i120X+52920X +I260OX +2100X +252X +21X +X 

2 3 4567 

6o4800X-(-12096ox +30240X ^-4200X +420X +28X -t-x 

2 3 4 5 h 7 

6o4800+42336oX+ 1J^1120X +294OOX +4200X +420X +28x +x 

2 3 4567 

166320OX+2268COX +65 520X +711OX *530X +35X +X 

2 3 IT 5 6 7" 

1663200+1058400X+31752CX +58800X +7520X +630X +35X +x 

2 3 4567 

399i6&3X+33 264ox -H3608OX +IOO8OX +924x +42X +X 
— 2 3 45 ~ T ~ 7 

3991630+2328J48OX+635040X -1-105840X -ill76ox -t882X -J-42X +X 

2 3 4-567 

864864ox+ 332640X •!-277200X -}-3.0080X +1512X +>t2X + X 

2 345 6^ 

86486403-465696ox+1i6424ox +1764oox +1764ox +U76x h49x +x 
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3 5 7 

G = 17297g30X+33 ^^0 0X -f- 30gHx +2X 

77 “ 2 3 5 1 

17297230436Ur;540X+i99584ox +2772C'OX -1-2320OX +1512X +56X +x 
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APPEITOIX 1(c) 



•-sT 

FADE' /J»PPxQXBIATION OF 1-e 

(Note: all listed complex roots liave complex conjugates which are 
not listed hut are represented by ’the open parenthesis,) 

X = sT 



G = 0 

00 



G 


X 


01 




G = 


X(X-2) 


02 


2 


G 


x(x-i.5-hii.956)() 


05 


6 


G = 


x(x- 2. 785 )(x-. 607 + 02 . 872 ) 0 


04 


24 


G 


X 


10 


X+1 


G = 


2X 


11 


X+2 


G 


x(x- 6 ) 


12 


2(X+5) 


G 


X(X-5+.15.875)() 


15 




G 


x (x- 5 . 458 ) (x-l . 281 +j 4 . 520 ) ( ) 


l4 


24(x+5) 


G 


X(X+2) 


20 


(X+l+ol)() 


G 


x(X4-6)_ ^ 


21 


(X+2+jl.4l4)() 


G 


12X 


22 


(X+5+jl. 752)0 


G = 


x(x~5+o7.i4i)() 


25 


Ix^4^j^ 


G 


X (X- 9 . 646 ) (X-l . 176 + o5 . 994 ) ( ) 


24 


“(X+5+j2.256Trr'^ 
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G = x( x-!-:ur>^ -.i i>936)( ) 

50 (x+l. 596 )'(x+. 702+0 1.807) d 

G = X(x+ 5 +j 5 o 875)() 

51 ■(X+27526) (X+l i 687-1- j 2 . 509 ) ( ) 

G = X( X+5+.1 7. l^l)() 

52 (x^ 37658 ) (x+2 . 681+0 5 . 050 ) ( ) 

G - 2X(X+. 1 7.746)() 

55 (X+C644)(x+5.678i-o5.509)() 

G = XiX--19 . 157 ) ( X- . 422+.i6 .080 ) ( ) 

54 (X+9.648) (X+4. 676+05. 915) 0 

G = X(X+2,785)(X+.607+.12.872)() 

40 (x+1 . 729+0 . 889 ) (x+ . 270+02.505 ) () 0 

G = X (x+5 . 4 58 ) (X+1 . 28l+.i 4 ♦ 520 ) ( ) 

41 l[x+2T7o4-fo 1 . 162 ) (x+l . 256+3 5 . ^57 ) ( ) ( ) 

G = X (X +9 . 648 ) (x-t-l . 176+0 5 .99U ) ( ) 

42 "(3^5 . 779+0 1.580) (X+2 . 221+3 4 a6o ) 0 0 

G = X (X+19 . 173 ) (X-f . 422 H-. 1 6 . 608 ) 0 

45 ^ . 787+ Ji. 567] (x+S". 215+3 ^ . 7751 0 0 

G • = to (x+.i 6.481)0 

44 (x+5 . 792+01 . 754 ) (X+4 . 208+ j 5 . 515 ) 0 0 

G = X( X-.i 84+.15.44 o)(x+ 2.6S4 +.1 1.704 ) 0 0 

50 Ix+ 2 .i 6 i) (X- .240+05.128) (x+1. 649+3 1.694) 00 

G = X ( X+ . 514 +o 4.844) (x-t-4 . @6-!-,l 2 . 95 2 ) ( ) ( ) 

51 7 x+ 5 . 257 ) (X+.705+J4.260 ) (x+2 . 678+0 2 . 181 ) ( ) (■) 

G = X(X+.409+.15.797)(X+7.0914-. 14.95‘ 0()0 

52 (x^ 47^2 ) (x+1 . 077+05.158 ) (X+5.690+J 2 . 571 ) 0 0 

G = X(X+.22 4+. i6.242) (x+9. 776+08. 755 ) ( ) ( ) 

55 ^5 . 277 ) (X+2. 665+0 5.884) (x+4 . 697+0 2 . 909 ) 0 0 

G = X f X+ . 0 58-} .1 6.527) (X+9 . 942+.1 16 . 698 ) () () 

54 “(x 737288 ) (x+5 . 656+36.544) (x+5 . 70l+j 5 . 210 ) ( ) () • 

G = 2x (x+,1 6 . 506) .(X +.U9... 5oo:):i:io 

55 (X+7.295) (X^J+.649+j7.l42) (x+6.7o4+j5.li85) ( ) ( ) 

G = X( X+5.555) (X-.874+.15.888) (X+2.097 +.1 2.892) O () 

60 Tx+2 . 562+0 . 858 ) (X- . 805+J 5.698) (X+1 . ) 4^2 +o 2 . 454 )()()() 
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G 

61 

G 

62 



x(x±s^gn(Xr^^±AiL 3 J!iKx±^ ^01 ) ( ) ( ) 

ix-o . ^ 24-i-j 1.047) (X-!- . 5 » 005 ) (X-!-^5 S+o 5 • 10270110 

X(X -!8.977) (X- .025-!-,15e8S7 ) ( X^4^.5564-■1 6 . 64 9)()() 
(X-{4.454-(-jl.2l7) (X+1.082i-j6.025) (X-^7r5B4+j5.639y( )(K) 



G = X(X -!-.048-!-.i6.i 86) (X-:-1'5.5 -56) (X-i- 5.284-i-.i9. 517) 0 0 

63 (x+5 . 471+0 1 736^ ) (x-^-2 . o6l-;- j 6 ^ 7 ) (X-t-4. 467-1- j 4.104) ()()() 

G = X (X-}20 . 598) (X-t-.026-f.l6. 2 75) fe-Jl.675-H 12. 831)0 Ci 

64 Tx-:'<7^!82-!-jl.499) (x-i-3.048-i-j7.652) (X-i-5 .469-i-j 4.520) ()()() 

G = X(X-^37.825) (X-i-.006-H.l6.287) (X-^2.082-i-.U4.770) 00 

65 (X-H7.491-J-0 1.622) (x-i-ii.039H-j8.346) (x-i-6.470-t-j4.900) ()()() 

G = 84x (X -Ki 6. 285 ) (XH-.114 .160 ) O O 

66 rxT87497H-jl . 735 ) (X-h5 . 032-i-j8 . 985 ) (x-i-7 . 47l-i-j 5 . 252 )()()( ) 



G = X(X-!-3. 6 o5-Kll.577HXH-1.4oo-i-.i3.700)(X“1.504-i-.i4.3 0 6)000 

70 1^2 . 759 )lx-iT4o6^i-j4 . 225 ) (XH-1 . l47-Hj 3.124) (x+2.380H-ji.6^) ( ) (TO 

G = X(x-!-5.684H-j2,443l (X-.8.30-O5.587 )(X-H2.l46-t -j5 

71 Ix77^B39lIx- .509-0 5T^5T(x^ (x+3.45iTj2.oi2)()()T) 

G = X(x~.245-o6.l 32lLXH:8.l42-!-j 3.631) (X-t-2.6o4H-j7..338) ()() ()^ 

72 Ix-tJT.BfS) (XH-.445-i-j6.837) (Xh -4 . 470 H-j 2.325) (Xh -3 . iK+jOlS ) ()()() 

G • = X(X--.032H-j6.2itG)(XH-2.862H-j9.6l4)(XH-ll.l70H-j5.438)()Il() 

73 90lT(X-Hl.4i8H-j7. 809) (X-!-5.486H-j2. 598) (XHii.l45+j5. 19^70170 

G = X(X-!-o012H-j6.287) ( X-!-2. 500-i-.ill.596)(X-n4.988H-.l8.625)O 00 

74 rx-H5 . 668 ) ^x-h 2 . 206-Hj8V732 ) (x-tit.928H-j4.726) (x-H7.552H-j4.567yO”OT) 

G = X(Xl^l±i^M)j:X;!Aj.5-19±ji2, 

75 (x?r^-0 3 . 071 ) ( x-t7 . 506 ) (xh-3 . 533h- j 9 . ii50 ) (x-i-6 . i42h- j 6 . i84 )()()( ) 



G 

76 



623)0()() 



G = X (X-O 6 . 285 ) (XH-.i 12 . 958 ) (XH-.i 36 . 13,9 ) ( ) O 0 ' 

77 W^^ 57 (X-i- 9 . 516-1-j 3. 478 ) (XH- 5 . 37 lH-jlO. 84 i) (x-i 8 .l 40 H-j 7 . 034) ()()() 
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JOB0608Ff SEILER* R.Wo 
PROGRAM FACTOR 

THIS PROGRAM FACTORS N POLYNOMIALS* WITH REAL COEFFICIENTS* 
OF ORDER 30 OR LESS, THE CALLING ARGUMENTS ARE 
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901 FORMAT (15) 

902 FORMAT ( IX* F9, 5 * 6F10, 5 ) 
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PROGRAM SPEC! 

THIS PROGRAM CALCULATES AND PLOTS THE SPECTRAL WEIGHTING 
FUNCTION OF I TRANSFER FUNCTIONS# CHARACTERIZED BY AN N-TH 
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